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STRUCTURE AND TAXONOMIC DIVERSITY OF BENTHIC 
DIATOM ASSEMBLAGE IN A POLLUTED MARINE ENVIRONMENT 
(BALAKLAVA BAY, BLACK SEA)‘ 


ALEXEI PETROV, ELENA NEVROVA, ANNA TERLETSKAYA, 
MIKHAIL MILYUKIN & VICTOR DEMCHENKO 


Abstract. The study examined changes in the structure and diversity of benthic diatom assemblages from 16 sampling stations 
covering ca 0.3 km? of the bottom of industrially polluted Balaklava Bay (Crimea, Black Sea). We analyzed data on diatom 
abundance and species richness along with measurements of 15 variables: depth, sediment size structure, TOC, metals (Zn, Ni, 
Cu, Pb, Hg, Cd, Cr, Mn) and organochlorine toxicants (PAHs, pesticides, PCBs). Based on cluster and MDS analyses, two groups 
of stations were distinguished at 30% similarity level. The first group comprised stations in the inner, most polluted part, and 
the second corresponded to the outer part of the bay. Diatom assemblages of a certain structure are formed within each group 
of stations. The species most responsible for similarity/dissimilarity within/between each of the distinguished complexes were 
determined and are listed. In total, 183 species representing 58 genera were identified. Species richness was highest in the genera 
Nitzschia Hassal (28 species), Amphora Ehrenberg (19), Navicula Bory (16), Cocconeis Ehrenberg (15), Diploneis Ehrenberg ex 
Cleve (11) and Lyrella Karayeva (6). Our analyses showed that the combination of depth, percentage of sandy fraction, and the 
content of TOC, pesticides, PCBs and mercury exerted the most significant impact on diatom assemblage structure and diversity. 
Changes in benthic diatom taxonomic diversity were evaluated for each station on the basis of average taxonomic distinctness 
(A*) and variation of the TaxD Index (A*). The taxonomic structure of diatom assemblages in Balaklava Bay is less diverse than 
the average level expected from the Crimean regional inventory (539 species). The A* value (80.89) was lowest for the most 
polluted part of the bay, where species-rich branches of taxonomic tree structure prevailed. The diatom assemblage in the part 
less impacted had a higher A* value (82.86), closer to the expected average mode (83.67) for the entire Crimean coast. Overall 
the results highlight problems with the use of taxonomic diversity measures to detect the impact of adverse environmental 
conditions based on benthic diatom assemblages. 
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INTRODUCTION 


Benthic diatoms (Bacillariophyta), the most abun- 
dant and ubiquitous group of the microphytob- 
enthos, have the highest population densities and 
species richness in coastal ecosystems. They are 
closely associated with certain biotopes directly 
subjected to environmental influences and sensi- 
tive to changes in the levels of various pollutants. 
These features recommend this group as a reliable 


* Dedicated to Dr. Kurt Krammer on the occasion of his 85" 
birthday 


indicator of eutrophication and anthropogenic pol- 
lution in comprehensive monitoring of aquatic eco- 
systems (Sladecek 1986; Sathya & Balakrishnan 
1988; Watanabe et al. 1988, 1990; Descy & Coste 
1991; Warwick 1993; Barinova et al. 2000, 2006; 
Gimenez-Casalduero 2001; De la Rey et al. 2004; 
Tavassi et al. 2004; Birkett & Gardiner 2005; Ni- 
kitina & Shkundina 2009). 

Most studies assessing biodiversity in the con- 
text of ecological monitoring and bioindication in 
marine and freshwater ecosystems deal with abun- 
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dance and species richness patterns. A variety of 
univariate numerical diversity indices combining 
species richness, abundance and evenness have 
been used to assess the effect of environmental 
degradation on biota, especially benthic assem- 
blages (see Magurran 2004; Karydis 2009 for 
review). Traditional diversity measures based on 
species richness and evenness (e.g., Shannon’s H’, 
Pielou J) often have disadvantages in assessment 
of biodiversity change on large spatial scales and 
long time scales. In the last decade or so a newer 
biodiversity measure emphasizing the average 
taxonomic relatedness between species in a com- 
munity have been developed and applied for bio- 
monitoring (Warwick & Clarke 1998, 2001). 

Measures based on the taxonomic structure 
of assemblages differ from more conventional 
diversity indices by incorporating the degree to 
which species are morphologically and thus evo- 
lutionarily related (according to the classical Lin- 
nean classification system for animals). Here we 
need to mention that the established taxonomic 
relatedness among diatom species is based pri- 
marily on frustule morphology and much less on 
genetic affinity. The premise is that an assemblage 
of closely related species must be regarded as less 
diverse than an assemblage of the same number 
of more distantly related species (e.g., in which 
the species belong to different taxonomic classes) 
(Warwick & Clarke 1998; Leonard er al. 2006; 
Leira et al. 2009). 

These biodiversity assessment indices have 
proved useful in several studies of different groups 
of organisms in different regions of the world (War- 
wick et al. 2002; Ellingsen er al. 2005; Leonard 
et al. 2006; Ceschia et al. 2007), but apart from 
a few studies focused on macroinvertebrate fauna 
(Heino et al. 2005, 2007; Campbell et al. 2007) 
the value of taxonomic distinctness measures as 
a means of expressing anthropogenic effects on 
the microphytobenthos is almost untested. 

Methods and indices for water quality assess- 
ment based on aquatic assemblages are rather 
widely applied in freshwater ecology (Stoermer 
1978; Tilman 1997; Kelly 1998; Dixit et al. 1999; 
Harding er al. 2005; Heino et al. 2005; Sonder- 
gaard et al. 2007), but there are very few studies 


relating to the use of marine benthic diatoms for 
bioindication of anthropogenically impaired bi- 
otopes (Dickman 1998; Izsak et al. 2002; Petrov 
et al. 2005, 2008). Most studies concerning Black 
Sea benthic diatoms deal mainly with taxonomy 
and floristic analysis, species composition and sea- 
sonal dynamics of diatoms (Guslyakov et al. 1992; 
Guslyakov & Nevrova 1998; Ryabushko 2008 for 
review). Fewer of them address the structure of 
the diatom assemblages and assess diversity fea- 
tures in relation to the influence of a broad set 
of environmental factors (Nevrova er al. 2003, 
2008; Ryabushko er al. 2003; Petrov et al. 2004, 
2005, 2007, 2008; Goldin 2009). Such compara- 
tive studies of diatom assemblage structure (es- 
pecially those considering the different levels 
of taxonomic resolution) can provide additional 
information about the influence of environmental 
changes and technogenic impacts on ecologically 
relevant measures of diatom diversity in Black 
Sea coastal regions. 

Our study assessed the variation of quantita- 
tive characteristics and taxonomic structure (ag- 
gregated to higher hierarchical levels) of benthic 
diatom assemblages resulting from environmental 
heterogeneity in Balaklava Bay, SW Crimea. We 
also analyzed differences in taxonomic distinct- 
ness indices between diatom assemblages along 
a gradient of permanent technogenic impacts (by 
several heavy metals and organochlorine com- 
pounds) in the bay bottom area. 


MATERIALS AND METHODS 


Data on benthic diatoms, including a list of species 
and their abundances in surface sediment samples to- 
gether with associated environmental information, and 
on bottom sediment chemistry, were obtained from 
the comprehensive ecological surveys carried out in 
Balaklava Bay (Black Sea, SW Crimea), in October 
2006. Sixteen stations representing the bay bottom at 
4-27 m depth were investigated (Fig. 1). Duplicate 
samples taken by tube (surface area 16 cm?) from the 
upper (1-4 cm) layer of soft sediments were collected 
at every station by a diver. 

Chemical analyses of inorganic and organic con- 
taminants in sediment samples included measurements 
of 14 substances: metals (Cu, Zn, Ni, Pb, Cr, Cd, Ag, 
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Mn, Hg), total PCBs, chlorinated pesticides and PAHs. 
The percentage content of silt+clay fractions and total 
organic carbon (TOC) of the sediments were also de- 
termined. 

The sampled sediments were prepared for chemical 
analysis according to the ISO 11464:2006 standard 
method. The samples were air-dried, sieved and homog- 
enized. Content of metals (besides Hg) in sediments was 
determined by graphite (MDL 0.005-0.05 mg/kg dry 
wt.) and flame (MDL 2.0-15.0 mg/kg dry wt.) atomic 
absorption spectrometry (AAS) following microwave 
digestion in Teflon pressure vessels with concentrated 
HNO;+HC1 (3:1). Total Hg in sediments was determined 
by cold vapor AAS. 

Organic contaminants (pesticides, PCBs) in sedi- 
ments were determined by GC/MS (MDL 1.0 ug/kg 
dry wt.) and capillary silica column GC/ECD (MDL 
0.05 mkg/kg) following Soxhlet extraction for 16 h with 
hexane/acetone (1:1). Total concentrations of PCB 
homologues as tetra-, penta-, hexa-, hepta-chlorobi- 
phenyls were determined. The extracts were cleaned 
by column chromatography and washing with H,SO, to 
remove interfering compounds. PAHs were determined 
by HPLC/UV (HP 1050/DAD, MDL 10-20 mkg/kg dry 
wt.) in reversed-phase mode. Standard solutions of indi- 
vidual substances or mixtures were used for calibration 
of equipment. Sediment grain-size composition as sandy, 
silty and clay fraction percentages was determined by 
wet sieving and gravimetric sedimentation. 

The benthic diatoms of all 16 sediment stations were 
analyzed. Duplicate diatom samples were examined for 
abundance, species richness and taxonomic structure. 

Only entire intact cells with visible living chloro- 
plasts were counted. All broken cells were rejected from 
microscope processing and were not used for quantita- 
tive and taxonomical analysis. Diatom species were 
identified to species or variety level (later referred to 
as species) using the taxonomic classification from 
Round er al. 1990, with additional information (Four- 
tanier & Kociolek 1999, 2003, 2007; Witkowski et al. 
2000). The diatom species were counted and expressed 
per 1 cm? of seabed. Species found only on permanent 
slides (i.e., rare or solitary species) were not included 
in the quantitative calculation of total density of dia- 
toms. These species, however, were included in the 
matrix (species density vs stations number) for further 
quantitative uni- and multivariate statistical analysis. 
In the matrix the density of such species was indicated 
as a conventional minimum value of 10 cells - cm”? 
(Petrov er al. 2005). 

Sample processing involved preliminary treatment in 
an ultrasonic bath for 20 min, followed by the standard 
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Fig. 1. Schematic map of sampling stations in Balaklava Bay 
in October 2006 (Black Sea, SW Crimea). 


technique using HCl and H,SO, with the addition of 
K,Cr,0,. Cleaned diatom valves were mounted using 
Eljashev mountant for light microscopy (Proshkina-La- 
vrenko 1974). Light microscopy (LM) employed a Carl 
Zeiss Axiostar microscope equipped with a PlanAPO 
100x objective. 

For this research we used the updated Crimean re- 
gional inventory (comprising 539 marine benthic dia- 
toms species) based on literature and our own materials 
(Nevrova & Petrov 2008). 

Multivariate statistical analysis was applied to eval- 
uate the spatial distribution of pollutants across the bay 
and to assess its influence upon the structure and taxo- 
nomic diversity of diatom assemblages. Data processing 
was done with the PRIMER v5 package (Carr 1997; 
Clarke & Gorley 2001). The environmental gradients 
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along the bay bottom were assessed from the results of 
PCA analysis. Cluster analysis and nMDS ordination 
of environmental factors were employed to distinguish 
the grouping of stations in relation to different levels 
of anthropopression (Clarke 1993; Carr 1997). The 
significance of differences between a priori separated 
groups of stations within the bay bottom area was tested 
with permutation/randomization ANOSIM (analysis of 
similarity). Similarity between stations was assessed by 
complete linkage of the Bray-Curtis similarity measure 
(for double square root-transformed biotic) and Eucli- 
dean distance (for non-transformed normalized abiotic) 
matrices (CLUSTER and MDS routines). The SIMPER 
routine was applied to reveal the most significant spe- 
cies (Clarke & Warwick 2001). Among such species the 
indicative ones are mainly responsible for the average 
similarity within the distinguished diatom assemblages, 
and discriminating species are mainly responsible for 
the differences between them corresponding to each of 
the considered locations. 

Spearman rank correlation coefficients (p) were used 
to detect the combination of environmental variables 
giving the best match for high similarities (low rank) in 
the biotic (abundance data) and abiotic matrices, that is, 
to recognize the set of abiotic factors best explaining the 
spatial differences in benthic diatom assemblages across 
the surveyed bottom area (BIOENV routine) (Clarke 
& Warwick 2001). 

Taxonomic distinctness indices (TaxDI) were used 
to estimate the diversity features of the diatom assem- 
blages of the bay (Clarke & Warwick 1998; Warwick 
& Clarke 1998). 

An initial binary presence/absence biotic matrix was 
constructed, and later the matrix data was aggregated 
into 7 taxonomic levels (from intraspecific to division). 
The taxonomic analysis was based on calculation of 
two taxonomic diversity indices. Average taxonomic 
distinctness AVTD (A') reflects the mean path length 
through the assemblage taxonomic tree (to the phylo- 
genetic common node) connecting every pair of spe- 
cies randomly selected from the list. The AVTD index 
characterizes the vertical taxonomic evenness along the 
phylogenetic tree of the given community. Variation in 
taxonomic distinctness VarTD (A*) describes the vari- 
ability (variance) of these pairwise path lengths (@ij) 
between each pair of species i and j, around their mean 
value A*. The VarTD index characterizes the horizontal 
unevenness of the taxonomic tree from the conventional 
averaged level. In other words, VarTD characterizes 
various representativeness of lower taxa (or successive 
ratios of lower taxa number to higher ones) for each 
of the branches composing the whole hierarchical tree 


Table 1. Environmental variables and variation ranges of their 
values used to characterize Balaklava Bay stations (samples 
were taken from upper 0—4 cm of bottom sediments). 
























































Environmental variable Variation range 
Depth, m 4.0-23.0 
Sandy fraction % 0-60.0 
Silttclay % 3.0-92.2 
TOC, % 1.3-6.6 
Organic toxicants, mkg : kg! dry weight 

PAHs 900.0-26000.0 
PCBs 2.0-435.0 
Pesticides 5.0-93.0 
Metals, mkg - kg! dry weight 

Cu 34.0-350.0 
Cd 0.1-0.6 

Zn 53.0-600.0 
Cr 2.5-67.5 
Mn 175.0-470.0 
Pb 51.0-1500.0 
Hg 0.3-2.0 

Ni 9.0-44.0 








(Clarke & Warwick 1999). The combination of AVTD 
and VarTD in a 2D scatter plot provides a statistically 
robust summary of the taxonomic relatedness patterns 
within the diatom assemblages. Calculated values of At 
and A* are usually located within the simulated 95% 
probability funnel converging towards increasing spe- 
cies number. The probability funnel (based on 1000-fold 
random pairwise selections) is created for a range of 
randomly drawn subsets (m = 20, 30, 40, ... 500, etc.) 
of species (or higher taxonomic categories) from the 
whole regional inventory, and a plot of the resulting 
upper and lower limits on a graph of A* or A* against 
m is constructed. The positions of real TaxD index values 
on such plots corresponded to certain stations, allowing 
comparison of deviations of the distinctness values (i.e., 
taxonomic structure) from the expected TaxD level cor- 
responding to the whole Crimean regional inventory. 
Points located within the probability funnel indicate 
that the taxonomic diversity of the corresponding sites 
falls within the expected range revealed from diversity 
analysis of the entire biogeographical region. Values 
outside the limits represent a biodiversity measure that 
departs significantly from what might be expected due 
to the presence of any pronounced (and long-term) en- 
vironmental perturbation (Warwick & Clarke 1998). 
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RESULTS AND DISCUSSION 


Measurements of the environmental parameters 
showed that the depth, grain size of bottom sedi- 
ment, TOC content and concentrations of toxic 
compounds vary considerably across the bay 
bottom (Table 1). There is a pollution gradient 
of heavy metals and organochlorine compound 
concentrations. 

Principal component analysis revealed a dis- 
tinct trend along the first three PC axes. PC1 (ac- 
counting for 32-36% of total explained variance) 
is associated with the gradient of heavy metals 
(Zn, Pb, Cu, Cr, Cd), while PC2 (27-28% of 
variance) is associated mainly with changes in 
TOC level and percentage of silty/clay fraction 
in the sediment. PC3 (12-14% of explained vari- 
ance) is associated with the distribution of PCBs 
and pesticides. 

Cluster analysis (based on similarity matrix of 
15 considered variables) divided the Balaklava Bay 
bottom into two well-distinguished zones, gener- 
ally corresponding to the inner (shallower) part and 
the outer (deeper) part and bay entrance (Fig. 2). 

Average values of abiotic variables recorded 
from the bottom sediments of those parts of the 
study area are given in Table 2. The content of 
most of the pollutants (except for Pb, PCBs and 
PAHs) was higher in the inner part ofthe bay. The 
most significant differences were between core 
stations from inner (station 8) and external (station 
16) parts of the bay. 

Diatom diversity changes were assessed versus 
the degree of anthropogenic impact on the basis 
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Fig. 2. Dendrogram representing the relative similarity of sta- 
tions (based on Euclidean distance measure for 15 key envi- 
ronmental variables). Dotted line indicates the integration level 
(ca 5.5) of sampling stations into two spatial groups: 1 — inner 
part of bay, 2 — outer part. 


of detailed data on sediment chemistry and other 
environmental parameters across Balaklava Bay. 

The updated Crimean inventory at present 
contains 539 benthic diatom species aggregated 
into 100 genera, 47 families and 25 orders. Micro- 
scopic analysis revealed the presence of 192 taxa 
belonging to 183 species, 58 genera, 38 families 
and 21 orders in Balaklava Bay (Table 3). Species 
richness was highest for genera Nitzschia Hassal 
(28 taxa), Amphora Ehrenberg (19), Navicula 
Bory (16), Cocconeis Ehrenberg (15), Diploneis 
Ehrenberg ex Cleve (11) and Zyrella Karayeva (6). 
Paraliales, Rhabdonematales, Thalassionematales 
and Toxariales were the most species-poor orders, 
with only one species recorded in each. The most 
abundant diatom species of Balaklava Bay are 
listed in Table 4. 

In addition to environmental gradient anal- 
ysis, spatial grouping of stations based on abun- 


Table 2. Mean values of environmental variables in bottom sediment corresponding to two parts of Balaklava Bay. Data for 
most polluted (station 8) and least polluted (station 16) core stations also indicated. 






























































Variables 
Location | 24 | Cu | Ni | Pb | Cr | Cd | Ag | Hg | Mn SI eg shen PAHs | TOC = 
Inner part | 206.5| 146.8| 33.1| 224.5] 34.0| 03 | 03 | 06 | 365.1| 842] 15.6| 37590] 23 | 809 
Outer part | 185.5] 130.7; 162] 317.6] 128] 0.1 | 03 | 1.1 1225.21 1227| 9.3|10470.0| 48 | 363 
St. 8 600.0| 350.0] 44.0 1500.0] 67.5] 05 | 04 | 08 | 325.0| 435.0/ 93.0] 27000 4.2 | 72.0 
St. 16 53.0; 338] 10.51 sı2| 13.8] 01 | 02 | 04 | 1920] 79.0] 12.0| 1200.0) 28 | 202 





Content of metals indicated in mg : kg! dry weight, organic toxicants in mkg - kg! dry weight of bottom sediment, TOC and 


sediment size fraction in percents. 
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Table 3. Composition of benthic diatom assemblages in Balaklava Bay. 



































































































































Taxon Sp. | Taxon Sp. 
THALASSIOSIRALES Gleser & Makarova 1986 RHOICOSPHENIACEAE Chen & Zhu 1983 

THALASSIOSIRACEAE Lebour 1931 Rhoicosphenia Grunow 1860 1 

Thalassiosira Cleve 1873 7 | CYMBELLACEAE Greville 1834 

STEPHANODISCACEAE Glezer & Makarova 1986 Cymbella Agardh 1830 2 

Cyclotella (Kützing) Brebisson 1838 4 | GOMPHONEMATACEAE Kutzing 1845 

Stephanodiscus Ehrenberg 1845 Gomphonema Ehrenberg 1832 2 

MELOSIRALES Crawford 1990 ACHNANTHALES Silva 1963 

MELOSIRACEAE Kützing 1844 ACHNANTHACEAE Kutzing 1845 

Melosira Agardh 1824 1 | Achnanthes Bory 1822 3 

HYALODISCACEAE Crawford 1991 Planothidium Round & Bukhtiyarova 1996 2 

Hyalodiscus Ehrenberg 1845 COCCONEIDACEAE Kutzing 1845 

Podosira Ehrenberg 1840 Cocconeis Ehrenberg 1837 15 
PARALIALES Crawford 1990 NAVICULALES Bessey 1907 

PARALIACEAE Crawford 1988 BERKELEYACEAE Mann 1990 

Paralia Heiberg 1863 1 | Parlibellus Cox 1988 3 

COSCINODISCALES Round & Crawford 1991 DIADESMIDACEAE Mann 1990 

COSCINODISCACEAE Kiitzing 1845 Luticola Mann 1990 1 

Coscinodiscus Ehrenberg 1839 3 | SELLAPHORACEAE Mereschkowsky 1902 

HEMIDISCACEAE Hendey 1937 em. Simonsen 1975 Fallacia Stickle & Mann 1990 3 

Actynocyclus Ehrenberg 1837 1 | PINNULARIACEAE Mann 1992 

TRICERATIALES Round & Crawford 1990 Pinnularia Ehrenberg 1843 2 

TRICERATIACEAE (Shiitt) Lemmermann 1899 Caloneis Cleve 1894 2 

Pleurosira (Meneghini) Trevisan 1848 1 | DIPLONEIDACEAE Mann 1991 

Triceratium Ehrenberg 1839 1 | Diploneis Ehrenberg ex Cleve 1894 11 
PLAGIOGRAMMACEAE De Toni 1890 NAVICULACEAE Kutzing 1845 

Dimeregramma Ralfs 1861 1 | Navicula Bory 1822 16 
Glyphodesmis Greville 1862 1 | Haslea Simonsen 1974 1 

FRAGILARIALES Silva 1964 Trachyneis Cleve 1894 1 

FRAGILARIACEAE Greville 1835 Dickieia Berkeley ex Kützing 1844 em. Mann 1994 1 

Fragilaria Lyngbye 1819 Astartiella Witkowski, Lange-Bertalot & Metzeltin 1998 | 1 

Synedra Ehrenberg 1830 PLEUROSIGMATACEAE Mereschkowsky 1904 

Opephora Petit 1888 Pleurosigma Smith 1852 2 

LICMOPHORALES Round 1991 Gyrosigma Hassal 1845 4 

LICMOPHORACEAE Kiitzing 1845 Rhoicosigma Grunow 1867 1 

Licmophora Agardh 1827 4 | PLAGIOTROPIDACEAE Mann 1990 

RHAPHONEIDACEAE Forti 1912 Plagiotropis Pfitzer 1871 1 

Delphineis Andrews 1977 1 | STAURONEIDACEAE Mann 1991 

ARDISSONEALES Round 1991 Stauronella Mereschkowsky 1901 1 

ARDISSONEACEAE Round 1990 THALASSIOPHYSALES Mann 1991 

Ardissonea De Notaris 1870 3 | CATENULACEAE Mereschkowsky 1903 

TOXARIALES Round 1990 Amphora Ehrenberg 1844 19 








TOXARIACEAE Round 1990 





BACILLARIALES Hendey 1939 
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Table 3. Continued. 






























































Taxon Sp. | Taxon Sp. 
Toxarium Bailey 1854 1 | BACILLARIACEAE Ehrenberg 1831 
THALASSIONEMATALES Round 1990 Bacillaria Gmelin 1791 1 
THALASSIONEMATACEAE Round 1990 Hantzschia Grunow 1877 4 
Thalassionema Grunow ex Hustedt 1932 1 | Nitzschia Hassal 1845 28 
RHABDONEMATALES Round & Mann 1990 RHOPALODIALES Mann 1992 

RHABDONEMATACEAE Round & Crawford 1990 RHOPALODIACEAE (Karsten) Topachevsky & Oksiyuk 1962 
Rhabdonema Kützing 1844 1 | Rhopalodia Müller 1895 2 
STRIATELLALES Round 1992 SURIRELLALES Mann 1991 

STRIATELLACEAE Kützing 1846 ENTOMONEIDACEAE Reimer 1976 

Grammatophora Ehrenberg 1840 3 | Entomoneis Ehrenberg 1845 2 
Hyalosira Kützing 1844 1 | AURICULACEAE Hendey 1965 

LYRELLALES Mann 1991 Auricula Castracane 1873 1 
LYRELLACEAE Mann 1991 SURIRELLACEAE Kutzing 1845 

Lyrella Karayeva 1978 6 | Petrodictyon Mann 1990 1 
Petroneis Stickle & Mann 1990 1 | Surirella Turpin 1828 2 
CYMBELLALES Mann 1990 Campylodiscus Ehrenberg ex Kützing 1844 3 





dance and species richness (biotic matrix) was 
performed. Such an analysis shows the spatial 
heterogeneity in diatom assemblage distribution 
through the bottom area, and it distinguished 
two well-separated groups of stations at the 30% 


similarity level (Fig. 3). Similarly to the results 
of environmental analysis, the first group (stations 
1-8 and 10) corresponds to the inner part, and the 
rest (9 and 11-16) are near the outer part and the 
mouth of the bay. 


Table 4. Leading complex of mass species of benthic diatoms in Balaklava Bay, their maximum and average abundance and 


occurrence. 




















: Abundance x 10° cells - cm? | Occurrence 
Species i 
Average Max % 

Nitzschia compressa (Bailey) Boyer 59.8 180.0 100 
Nitzschia sigma (Kiitzing) W. Smith 30.2 90.0 100 
Thalassionema nitzschioides (Grunow) Mereschkowsky 6.0 33.8 100 
Cocconeis scutellum Ehrenberg 2:5 11.3 100 
Bacillaria paxillifera (O. Müller) Hendey 16.2 124.0 93.8 
Nitzschia coarctata Grunow 1.8 45.0 93.8 
Diploneis smithii (Brebisson) Cleve 6.4 22.5 93.8 
Grammatophora marina (Lyngbye) Kützing 4.9 22.5 93.8 
Nitzschia panduriformis Gregory 2:5 16.9 93.8 
Synedra tabulata (Agardh) Kützing 6.0 39.4 87.5 
Lyrella abrupta (Gregory) Mann 6.4 16.9 87.5 
Cymbella angusta (Gregory) Guslyakov 0.4 5.6 81.3 
Fallacia forcipata (Greville) Stickle & Mann 3.2 45.0 62.5 
Amphora coffeaeformis (Agardh) Kützing 2.8 33.8 62.5 
Caloneis liber (W. Smith) Cleve 5.6 33.8 53.0 
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SIMPER analysis revealed the most significant 
indicative and discriminating diatom species of 
each of the assemblages considered, presented in 
Tables 5 and 6. These results distinguished diatom 
assemblages of particular structure within each 
separated part of the bay. In the inner (polluted 
part), only 2 of 145 top-ranked species, Nitzschia 
compressa (Bailey) Boyer and Nitzschia sigma 
(Kützing) W. Smith, constituted ca 82% of total 
input into the average similarity within the assem- 
blages. In our interpretation these are the most in- 
dicative species of the assemblages inhabiting the 
highly polluted part ofthe studied bottom area. The 
diatom assemblage of the less polluted outer part of 
the bay included 5 top-ranked indicative forms (of 
the total 136): Grammatophora oceanica Ehren- 
berg, Bacillaria paxillifera (O. Müller) Hendey, 
Nitzschia compressa (Bailey) Boyer, Grammato- 
phora marina (Lyngbye) Kützing and Thalas- 
sionema nitzschioides (Grunow) Mereschkowsky. 
These species together constituted ca 50% of the 
total contribution to average similarity within this 
assemblage. However, average similarity between 
stations within the distinguished pollution-related 
groups was rather low, especially for group 2 
(22.3%). This suggests that environmental condi- 
tions (pollution level) did not significantly influ- 
ence the existing differences in the distribution of 
dominant species within each of pollution-related 
parts of Balaklava Bay. 

On the other hand, comparative analysis of the 
quantitative development of top-ranked discrimi- 
nating species in the two parts of the bay revealed 
high dissimilarity (83.4%) between assemblages. 
This indicates pronounced differences in key spe- 
cies abundance between the compared assemblages 
living under different levels of pollution. 

The division of all stations into two spatial 
groups according to levels of environmental fac- 
tors and information about the set of species cor- 
responding to each group can also be used to find 
the combination of variables that best matches 
the high similarity between the biotic and abiotic 
matrices. It can identify the set of abiotic factors 
best explaining the spatial differences in benthic 
diatom assemblage patterns across the surveyed 
bottom area. 








Similarity 














Fig. 3. Results of clustering (A) and nMDS ordination (B) 
analysis: grouping of stations from Balaklava Bay based on 
diatom abundance. 


Evaluation of the Spearman rank correlation 
showed that the highest value of the coefficient 
(Pmax = 9.612) corresponded to a combination of 
depth, percent of sand fraction, PCBs, TOC, pesti- 
cides and Hg. This set of factors gave the maximal 
matching with the difference in biotic parameters 
and presumably is the cause of changes in the 
structure and diversity features of the diatom as- 
semblage along the pollution gradient in the bay. 

Earlier work on the distribution of benthic dia- 
toms along environmental gradients in similarly 
polluted Sevastopol Bay (Black Sea, Crimea) gave 
a combination of seven key abiotic factors: depth, 
percentage content of sand fraction, PCBs, Eh, Cd, 
Mn and Zn (Petrov et al. 2005). The first three are 
the same in the sets of key factors influencing the 
abundance and diversity of diatoms in the Sevas- 
topol and Balaklava regions. The rest of the key 
variables for Sevastopol Bay (Cd, Mn, Zn, redox 
potential Eh) differ from set of key variables de- 
tected in the current analysis for Balaklava Bay. 
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Table 5. Contribution of the most significant indicative species in average similarity within benthic diatom assemblages in two 


compared locations of Balaklava Bay. 


























Species in compared locations N, cells - cm? S; Si(%) 
1 — inner part (145 sp.) — average similarity 45.7% 
Nitzschia compressa 9136.7 22.1 48.5 
Nitzschia sigma 4688.9 15.3 33.5 
Other species 18.0 
2 — outer part (136 sp.) — average similarity 22.3% 
Grammatophora oceanica 2092.8 4.3 19.3 
Bacillaria paxillifer 3381.8 2.6 11.8 
Nitzschia compressa 1932.3 1.9 8.8 
Grammatophora marina 1049.0 1.8 8.1 
Thalassionema nitzschioides 1210.4 1.2 5.4 
Other species 46.6 





Note: N, cells - cm”? — average abundance of i-th species in assemblage, S — similarity function, S;— absolute and S; (%) — relative 
contribution of i-th species in average similarity within the diatom taxocenotic complexes. 


Dickmann (1998) noted a significant correlation 
(P < 0.05) between high levels of sediment toxicity 
from metals (Cd, Cr, Cu, Ni, Zn) and decreasing 
species diversity and deformity of cell morphology 
of diatoms of the genera Navicula, Achnanthes 
and Fragilaria. 

The next part of this study used AvTD (A') and 
VarTD (A*) indices to compare possible deviations 
in diatom assemblage diversity across the bottom 
of the bay. We constructed 95% probability 2-D 
funnel plots for the simulated AVTD and VarTD 


values, based on the inventory of the Crimean ben- 
thic diatom flora, and points for each station were 
superimposed on the simulation funnel. The AVTD 
value (81.59) calculated for the whole Balak- 
lava Bay was below the expected average cor- 
responding to the entire Crimean region (83.67). 
The AvTD value corresponding to the inner part 
of the bay was lowest (80.89). The phylogenetic 
structure of the diatom assemblage in permanently 
polluted Balaklava Bay is shown to be character- 
ized by depressed taxonomic diversity, especially 


Table 6. Contribution of the most significant discriminating species in average dissimilarity between benthic diatom assemblages 


(83.4%) in two compared locations of Balaklava Bay. 



































Species in compared locations ; N, Selle seng” D D; (%) 
1 — inner part 2 — outer part 
Nitzschia compressa 9136.7 1932.3 15.9 19.1 
Nitzschia sigma 4688.9 888.0 8.9 10.7 
Bacillaria paxillifer 254.7 3381.9 5.3 6.4 
Grammatophora oceanica 67.0 2092.9 4.4 5.3 
Synedra tabulata 320.0 968.0 2:3 2:7 
Grammatophora marina 70.3 1049.0 2.2 2.6 
Thalassionema nitzschioides 132.9 1210.4 2.2 2.6 
Coscinodiscus sp. 1 0.0 724.3 1.9 2.3 
Other species 48.3% 

















Note: N, cells - cm”? — average abundance of i-th species, D — similarity function, D;— absolute and D, (%) — relative contribution 
of i-th species in average dissimilarity (83.4 %) between the diatom complexes. 
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Fig. 4. 95% probability funnel plots for simulated AvTD, 
At (A) and VarTD, A* (B) values, based on the total list of 
Crimean benthic diatom flora (502 species). Points indicate 
the observed A* and A* values corresponding to diatom as- 
semblages from each of 16 stations (1st group — stations from 
inner part, 2nd group — stations from outer part of bay). Dotted 
line indicates expected average values of TDIs for the whole 
Crimean region. 


at higher hierarchical levels (genus, family). The 
group of stations corresponding to the outer, less 
polluted part of the bay has a slightly higher At 
value (82.86), close to the expected average mode 
from diversity analysis of the Crimean inventory. 
Figure 4A is a funnel plot with the points indicated 
for the observed AvTD values corresponding to 
each of the 16 sampling stations. Most of the points 
are near or below the expected level but within the 
lower limit of the probability funnel. 

The lower positions of the station points on the 
AvTD plot (according to expected average mode) 
is evidence that diatom diversity structure at most 
of the sampling sites (except stations 13 and 15) 
is below expectations. Such diversity features are 
related to decreased taxonomic diversity (in spite 
of considerable species richness) and with minimal 
evenness of the whole taxonomic tree due to the 
rather high ratio of species-rich branches. The 
most topologically flat (or less vertically even) 
hierarchical structure of the taxocene at stations 
1, 2, 5 and 14 also signify the lowest taxonomic 


diversity. This result is likely due to elimination of 
some species-poor and mono-specific branches of 
the taxonomic tree under permanent environmental 
disturbances, including anthropogenic pollution 
(Warwick & Clarke 2001; Warwick er al. 2002; 
Leonard er al. 2006; Leira et al. 2009). 

In the same type of comparative analysis per- 
formed with the other taxonomic diversity index 
(VarTD), only the points of stations 2, 6 and 15 
were close to the expected mode of the corre- 
sponding funnel (A* = 308) (Fig. 4B). These re- 
sults mean that the aggregation pattern of initial 
taxa into higher hierarchic levels at these stations 
is similar to the expected average value of VarTD 
for the entire Crimean region. 

The taxonomic diversity of the rest of the sta- 
tions (except station 13) was characterized by wide 
variability of phylogenetic tree topology: the cal- 
culated A* values mapped on the plot are not only 
above the expected average level but also beyond 
the upper limit of the 95% probability funnel (see 
Fig. 4B). This suggests that the diatom assem- 
blages are characterized by uneven phylogenetic 
structure (i.e., both species-poor and species-rich 
arms in the taxonomic tree). Taxonomic structure 
can become asymmetric under the influence of 
different environmental factors changing consid- 
erably across temporal and spatial scales. In such 
cases the assemblage is usually characterized by 
pronounced taxonomic heterogeneity, becoming 
apparent at higher hierarchical levels. Note that 
the interposition of points on plots A‘ and A* corre- 
sponding to the same station usually have a mirror 
pattern (e.g., see points of stations 1, 5 and 13), 
where the expected average line can be the axis 
of symmetry (see Fig. 4A & B). 

Results from a similar analysis of taxonomic 
structure considering higher hierarchical levels 
(genus, family) are given in Figure 5. The posi- 
tions of sampling points on the probability funnel 
plots for three hierarchical levels of diatom di- 
versity in Balaklava Bay indicate a certain shift 
of points from symmetrical positions around the 
expected average (dotted line) corresponding to 
species level. The ellipses on the plots highlight 
such a shift of a sampling point cloud from the 
initial rather symmetrical position. 
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Fig. 5. Simulated AvTaxD (A*) values: probability funnel plots 
for three hierarchical levels of diatom diversity in Balaklava 
Bay. A — species level, B — genus level, C — family level. 
Ellipses indicate shift of sampling points from the initial sym- 
metrical position around the expected average (dotted line), 
corresponding to species level. 


A comparison of these data with the above re- 
sults for species level provides additional informa- 
tion regarding alteration of the vertical evenness 
of the diatom assemblage phylogenetic tree, that 
is, degradation of taxonomic structure under long- 
term environmental impacts in the bay. 

Univariate relationships between key environ- 
mental variables distinguished after Spearman co- 
efficient analysis (rank comparison of biotic and 
environmental matrices), and changes in the AVTD 
index, are presented in Figure 6. The curve of 
changes of A* values along the axes of environ- 
mental variables (TOC, sandy fraction) and pollut- 
ants (PCBs, pesticides, Hg) is bell-shaped. These 
newly found relationships show deviation of taxo- 


nomic evenness (AvTD index) from the average 
level, due to the influence of various environmental 
factors. Taxonomic diversity also increases versus 
initial levels due to the relative prevalence of 
species-poor branches in the community structure 
hierarchy. Such changes could be induced by in- 
creasing environmental disturbances in the biotope 
(Leira et al. 2009). However, further increases of 
environmental disturbances (including pollution) 
can lead to degradation of community structure, 
the disappearance of species-poor branches, and 
their gradual replacement by the species-rich arms 
in the taxonomic structure of the assemblage. 
Other authors have also noted that disturbances 
(natural or anthropogenic) of a biotope tend to 
eliminate species-poor arms, leaving those that 
belong to comparatively species-rich taxonomic 
groups (Clarke & Warwick 1998, 2001). Generally 
the taxonomic distinctness of degraded locations is 
usually significantly lower than that of relatively 
pristine locations for different groups of organisms 
(e.g., benthic nematodes, coastal fishes, echino- 
derms) (Warwick et al. 2002; Ellingsen er al. 2005; 
Leonard et al. 2006; Ceschia et al. 2007). 

The AvTD and VarTD indices are useful tools for 
integrated assessment of changes in diversity under 
environmental degradation (Izsak & Price 2001). 
They could be important components of monitoring 
schemes for marine coastal ecosystems, in which 
evaluation of benthic diatom diversity is considered 
as a key bioindicator (Leonard et al. 2006). 

Our study is one of the first attempts in world 
diatomology to make a comparative assessment of 
taxonomic diversity features of marine benthic di- 
atom assemblages under different levels of techno- 
genic pollutants. Our results revealed relationships 
between the values of TaxD indices and levels of 
key pollutants (heavy metals, PCBs, PAHs) in Black 
Sea bay sediments. An approach based on compara- 
tive evaluation of marine benthic diatom diversity 
can be recommended for practical monitoring of 
marine coastal ecosystems subject to permanent an- 
thropogenic impacts, along with other bioindication 
methods using the different groups of the freshwater 
microphytobenthos (Schoeman & Haworth 1986; 
Stevenson & Pan 1999; Harding et al. 2005; Heino 
et al. 2005; Karydis 2009). 
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Fig. 6. Relationship between key environmental variables and AvVTD index (A* values) for diatom assemblages from Balaklava 
Bay stations. A — depth, m; B — TOC, %; C — sand, %; D — PCBs, mkg - kg; E — pesticides, mg - kg; F- Hg, mg : kg. 


CONCLUSION 


Cluster and MDS algorithms applied for united 
analyses of 15 abiotic variables and diatom abun- 
dance matrices revealed two spatially separated 
groups of stations in Balaklava Bay: one corre- 
sponding to the inner part of the bay and the other 
near the outer part and the mouth of the bay. The 
content of the majority of the pollutants, including 
heavy metals and organochlorine toxicants, was 
significantly higher in the sediments of the inner 
part of the bay. 

Pronounced differences in the abundance, spe- 
cies richness and taxonomic diversity patterns of 
benthic diatom assemblages were correlated with 
the level of anthropogenic impacts in Balaklava 
Bay. Certain diatom assemblages corresponding to 
each of the two separated groups of stations were 
distinguished. The indicative and discriminating 
species responsible for similarity within each of 


assemblages and dissimilarity between them were 
distinguished and are listed. 

Evaluation of Spearman rank correlation coef- 
ficients (Pmax = 0.612) showed that a combination 
of variables — depth, ratio of sand fraction, content 
of TOC, pesticides, PCBs and Hg — was best cor- 
related with changes in biotic patterns, suggesting 
their impact on the structure and diversity features 
of the diatom assemblages. 

Based on assessment of TaxDI, the phyloge- 
netic structure of diatom assemblages in the per- 
manently polluted Balaklava Bay is generally 
characterized by depressed taxonomic diversity 
as compared with the expected average level cor- 
responding to the entire Crimean region. The value 
of the AvTD index was lowest for the inner part 
of the bay, for which the taxonomic tree of the 
diatom assemblage is characterized by the relative 
prevalence of species-rich branches. The struc- 
ture of assemblages from the less polluted external 
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part of the bay gave higher AvTD values close 
to the expected average mode corresponding to 
the Crimean inventory. The results suggest that 
habitat variability and the gradients of physico- 
chemical factors through the Balaklava Bay alter 
the abundance and erode the taxonomic structure 
of the diatom assemblage in conditions of long- 
term pollution. 

Taxonomic distinctness indices are important 
tools for further development of marine coastal 
ecosystem monitoring, in which benthic diatom 
diversity is a key criterion among the bioindica- 
tion approaches. 
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